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Abstract The main goal of this work is investigation 
of NADE in the cyclic universe scenario. Since, cyclic 
universe is explained by a phantom phase (oj < —1), it 
is shown when there is no interaction between matter 
and dark energy, ADE and NADE do not produce a 
phantom phase, then can not describe cyclic universe. 
Therefore, we study interacting models of ADE and 
NADE in the modified Friedmann equation. We find 
out that, in the high energy regime, which it is a nec- 
essary part of cyclic universe evolution, only NADE 
can describe this phantom phase era for cyclic universe. 
Considering deceleration parameter tells us that the 
universe has a deceleration phase after an acceleration 
phase, and NADE is able to produce a cyclic universe. 
Also it is found valuable to study generalized second law 
of thermodynamics. Since the loop quantum correction 
is taken account in high energy regime, it may not be 
suitable to use standard treatment of ther modynam- 



ics, s o we turn our attention to the result of (|Li et al. 
20091) . which the authors have studied thermodynam- 
ics in loop quantum gravity, and we show that which 
condition can satisfy generalized second law of thermo- 
dynamics. 
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1 Introduction 



Cosmological and astronomical obser vations such as su 



pernoyae type la observat i onal data (IPerlmuter et al 



1997t iReiss et al. I Il998t Ide Bernardis. et al. I 120001) 



and Wilkonson Microwave Anisotropic Probe (WMAP) 



(Astier et al. 12006 


JPerlmutter et al. Il999l 


Peiris et al. | 


20031 Riess et al. 12007) imolv that the 


universe is 



undergoing a period of accelerated expansion. Since 
normal matter can not give rise to accelerated expan- 
sion of the universe, Scientists came up with a solution 
which expresses that this expansion is a result of an am- 



biguous fluid name d Dark Energy ([Bridle et al. 112003 



Koght et al. 1 12003) 



The theoretical and experimental analysis suggest that 
the universe consist of %73 dark energy, %23 cold 
dark matter (CDM), and remnant matter is baryons 



( Fuiii and Maeda 1 120031 : ISpergel et al. 1 120031 ). Unfor 



tunately the nature and origin of dark energy are am- 
biguous up to now, but people have proposed some 
candidates to describe dark energy. Amongst the var- 
ious candidates of dark energy to describe accelerated 
expansion of the universe, cosmological constant (vac- 
uum energy), A, with equation of state (EoS) ui = — 1 
is located in central position. However, as it is well 
known, the cosmological constant proposal has two 
famous problems, fine-t uning problem and the cos 
mic coincidence problem ( Weinberg! 1 1989j Isteinhardt 
19971 ICarroll I l200ll Peebles and Ratra 1 120031 ). Some 



other of dark energy models suggest that dark en- 
ergy component can treat as scalar field with dy- 
namical EoS. In this scenario the evolution of the 
field is very slow, so that kinetic energy density is 
less than the potential energy density, and this give 



celeration (Frieman et al. 


1995 




Turner and White 


1997: 


Caldwell et al. 1998 


). Some of scalar field 



models are as chameleon field ([Khourv et al. 1 12004 ; 
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Mota and Barr ow 200 4J) , quint essence (Q-field) (jLiddle. 
19991 IPeebles and Ratra Ill999h . and phantom. In the 



phantom field scenario, the parameter of EoS is as 
us < — 1, due to existence of a negative kinetic en- 
ergy density of scalar field. It is well known that the 
phantom dark energy model suffers from two kind of 
problems, "Big Bang" singularity and "Big Rip" sin- 
gularity, where big bang is related to initial epoch of 
universe and big rip is related to a finite future sin- 
gularity. Since the space-time singularities are invidi- 
ous for theorists, models that avoid these singularities 
are attractive. One of these models na mes cyclic uni 



etlEhcfrist and Sloth 112004 iDutta et al. Il2010l) . Another 
attractive model to describe the nature of dark energy, 
within the framework of a fundamental theory origi- 
nating from some considerations of the feature of quan- 
tum gravity t heory, is c alled agcgraphic dark energy 
(ADE) model |Cai Il2007t) . The ADE assumes that the 
dark energy comes from the universe components fluc- 
tuation such as space-time and matter fluc tuation (for 



verse which have received huge attention ( Brown et al . 
20081 iBaum and Frampton ||2007t iBrown et al. II2008I ). 
Presence of p 2 term with a negative sign in Friedmann 
equation (which is used at studying cyclic universe) is 
an effective way to eliminate these singularities. In the 
cyclic universe scenario, where is based on the phan- 
tom dark energy model, the universe oscillates through 
a series of expansion and contraction. Universe in this 
scenario has a very high energy density at beginning an 
ending of the expansion, s o quantum gravity can not 



be ignored in these stages ( Ashtekar. 2006bl la). This 



evolution can be result from the modified Friedmann 
equation in the loop quantum cosmology (LQC). In 
LQC, the Friedmann equation has been modified to 



H 2 = 



P 



3m 2 



Pc 



(1) 



where H is the Hubble parameter, m 2 is the reduced 
planck mass (m 2 — — ^ = 2.44 x 10 18 GeV), p is the to- 
tal of energy density, p c is the critical energy density as 
p c = 4V37- 3 ™* = 0.82p p , where p p = 2.22 x 10 76 GeV, 
and 7 is the dimensionl ess Barbero-Immirizi parameter 



( Ashtekar et al. Il2006d) . We notice, this correction can 
solve the singularity problems as follow, when the to- 
tal energy density reaches the critical density, the uni- 
verse reaches the maximum at the end of expansion 
that is called "turnaround point", and universe arrives 



bounce there (Bovle. et al. 


~~ ~ "i " — ^ " w 

2004; Khourv et al. 


2001 


Steinhardt and Turok 2002 


; Bovle. et al. 2006) 


. We 



emphasize the ide a of cyclic universe was first intro 



duced by Tolman (jTolman Ill93ll 119341) 



An interesting attempt for probing the nature of dark 
energy, in the framework of quantum gravity, is the 
holographic dark energy (HDE). In the HDE model, 
dark energy is a dynamical evolving vacuum energy 
density th at can satisfy the p hantom behavior. Au- 
thors of ( Zhang et al. I |2007| ) have investigated the 



cyclic universe by HDE (and some in t eresting work 



about HDE have been done (Li 



2007aHbl ; IKarami and Fehri 1 12010 ; 



2004; 



Setare et al. 



Cohen et al. 



1999; 



further discussion we refer the reader to (Kim et al 



2008: iNozri and Azizi 1 120091; iKim Il2008t IKarami et al 
2008tlSetarell2010l) . In this model, the age of 



: universe is 

taken as the length measure instead of the horizon dis- 
tance, therefore the causality problem which appears 
in the HDE model can be avoided. The ADE model 
suffers from the difficulty to de scribe the matter dom- 
inant epoch. The authors of (|Wei et al. 1 120081 ) have 
introduced a new mechanism to overcome that prob- 
lem, which it is called new agegraphic dark energy 
(NADE) model, and its energy density is defined by 
Pa = 3n 2 m 2 ri~ 2 , where 3n 2 is introduced to parame- 
terized some uncertainties and r\ is conformal time and 
can be written as 



?? = 



da 



(2) 



where a is scale factor and H is well-known as Hubble 
parameter. Cyclic model of universe, due to avoiding 
singularity, and NADE, due to estimating a good ap- 
proximate of dark energy value and solving causality 
problem of HDE, have received hu ge interest. In pre- 
vious works, like ( Wei et al. 20081 ). it was only men- 
tioned that interaction NADE can produce phantom, 
however it was not explained that if NADE can pro- 
duce cyclic universe and how the quantities like dark 
energy density parameter, equation of state parameter, 
and deceleration parameter behave in cyclic universe 
by taking NADE as component of dark energy. In this 
work we motivated to take NADE as dark energy and 
investigate universe evolution in cyclic model. If it is 
able to stand in phantom area, it may produce cyclic 
universe, however it is not all of story. Although, pres- 
ence of this type of dark energy displays an accelerated 
universe, but the most important thing in cyclic uni- 
verse is that, after a while, the universe should enter in 
a deceleration area and recontract in turnaround point. 



2 High and low energy 

In this step, we focus on high energy regime. Certainly 
the quadratic term of energy density in the Friedmann 
equation can not be ignored, this term can play a very 



3 



impressive role in the evolution of universe. The modi- 
fied Fricdmann equation is as 



Assumption of interaction between p m and p\ may 
solve above problem. By including interaction, the con- 
servation equations are as 



3H Z = p 1 



Pc 



(3) 



where p is a combination of dark energy density p\, 
and matter density, p m . For dark energy dominant in 
high energy regime, the above relation can be rewritten 
as 



p m + 3H(1 +oj m )p m 
pA + 3H(l + uj A )p A 



-Q, 



(8) 
(9) 



3H 2 



PA 



1-^ 



where Q indicates interaction. Q is taken as Q = TpA, 
with r > 0, w hich means the r e is tr ansfer of energy 

(l2010h . We take Q as 
= Setting ADE as 

PA ° 



(4) 



so, according to this relation the parameter of dark en- 
ergy may be estimated as 



from pa to p m IShevkhi et al. 
Q = 3b 2 H(l + r)p A , where r 
the component of dark energy gives us an uja as 



(10) 



PA 

3H 2 



Pc 



J7 A is always larger than one, and it becomes very large 



PA 



(5) 



'A 

at turnaround point, so for having uja < 
2^/n^ 



T, we should 



have b (1 



r) > 



3n 



. This predicts a large value for 



tained value of another papers such as (jShevkhi et al. 
2010h . 



Q A is always larger than one, and it can be very large coupling constant b, while it is in cont rast to the ob 
when the universe approaches to the turnaround point. 
However note that, since we have p m , although it is 
very small, f2 A can not be infinite. This small value of 
p m does not allow p A to reach the exact value of p c , 
then the presence of p m can prevent the infinite value 

of n A . 

If we assume there is no interaction between these two 
components of the universe, conservation equation can 
be written as 



Now, NADE is taken as p A which that is the main case 
of this work. Because we use conformal time r\ instead 
of T, scale factor appear in relation (fTU)) , namely 



(11) 



uj A = -l + —dnl-b 2 (l + r). 
Sna v 



p m + 3i?(l + U) m )p m = 

p A + 3H{l + uj A )p A = 0. 



To have uj a < — 1, the coupling constant b should obey 
following relation 



(6) 



By taking ADE, which is defined as p A — 3n 2 T~ 2 where 
T denotes time, as dark energy component of universe, 
lo a is obtained as 



b 2 



> 



Sna 



(12) 



3// 



(7) 



Also, if we replace NADE instead of ADE, we arrive 
at same result but only scale factor, a, should be added 
in the denominator of second term on the right hand 
of relation ([7]). Since the latest term, in both case, is 
always positive wa never can be smaller than — 1. Since 
phantom type of dark energy is requisite for cyclic uni- 
verse (phantom energy density is getting larger with in- 
creasing scale factor so the total energy density reaches 
the critical energy density in turnaround point and uni- 
verse begins contraction), selecting ADE and NADE, 
without interaction, as a com ponent of dark en ergy in 
cyclic universe is not suitable ( Wei et al. 20081) . 



here, r has been ignored because of p rn -C Pa- If ^ A be 
in order of a 2 , we may obtain a convenient value for b 
in order it could give us phantom dark energy in this 
regime. 

In Fig.l wa parameter has been plotted versus r = 
^j- , for three different value of interaction coupling con- 
stant. In expansion phase, by passing time, the quan- 
tity r decreases and makes lu a larger. Since it still 
stands in phantom range, larger value of uj a shows that 
dark energy density grows up slower; in contrast to the 
contraction phase. Also, having stronger interaction 
coupling causes smaller value for uj a . 

If f2 A be in order of a^, where £ < 2, equation ([12)) 
can be valid in good approximation. Now suppose f2 A 
is in order of a 2 , for instance f2 A = aa 2 . From cqua- 
tio n (fT2l) we have b 2 > and by taking the value 

of |Shevkhi et al. Il2010l) for b, namely b 2 = 0.25, we 
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tion, which governs the NADE evolution of universe in 
high energy regime, can be attained as 





b 2 =0.5 




ir = l 

fr=2 



T - , - r- - i_ 



T5.0B ■ 



Fig. 1 equation of state parameter of dark energy has been 
plotted versus the dimensionless quantity r — , for three 
different value of interaction coupling constant. 



obtain the value of a as yfa < |n. 

Now, we want to obtain differential equation for f2 A . 

the NADE, p A is given as 



In 







-2fi 



(n 



36 2 (l + r)(-ft A 



(17) 



with attention to (|12p . it is clearly shown that f2 A = 
■g-fi A is positive. So one can realize that when the uni- 
verse is growing up (decreasing) f2 A is positive (nega- 
tive) which it displays that dark energy density parame- 
ter is increasing (decreasing) by passing time. By pass- 
ing the universe from turnaround point, f2 A decreases 
and puts the universe in low energy regime. 

Another useful cosmological parameter is decelera- 
tion parameter which can tell us if the expansion of 
universe is accelerating or not. Deceleration parameter 
is given by 



PA 



3^ 

T] 2 



(13) 



From the definition of dark energy density parameter, 



namely Q A 
acquired as 



The differential equation for ftj 



is 



-1 - 



H 



(18) 



substituting -Mg term in the above relation, we arrive 
at 



-2Q 




(14) 



where prime denote derivative with respect to M = 
In a. Taking the time derivative of modified Friedmann 
equation, and substituting that in the above equation, 
we obtain 



q = -i + |((i + wm)n^ + (i + WA)ni) 

-l + !(l+ WA )Q A (l-^). (19) 
Deceleration parameter has been plotted versus the 



H 

IP 



((l + a; m )fi!L + (l + WA)fii) ( ] 



2 J1 

Pc 



(15) 



since, in the high energy regime, fijj^ can be ignored 
against to A , therefore one can estimate 



(16) 



(for driving above equation we have used — = - 



x Pc an" 

J 1 1 1 k 

see ref . (jShevkhi et al. 1 120101 ) . The differential equa 





;r=0.5 




Zr=l 




b-=2 




Fig. 2 (expansion epoch) Deceleration parameter, q, has 
been plotted versus the dimensionless quantity 9 = for 
three different value of interaction coupling constant. 
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— ir =0.5 



By defining the energy density parameters in the form 
of 



i 2 =l 
Z) 2 =2 




Fig. 3 (contraction epoch) Deceleration parameter, q, has 
been plotted versus the dimensionless quantity x = 1 — 9, 
for three different value of interaction coupling constant. 



dimensionless quantity 9 = ^ for expansion and con- 
traction epoch of universe evolution in Fig. 2. and Fig. 3. 
respectively (with attention to this fact that uj\ stands 
in phantom range). In the plots three different val- 
ues for b are taken as b 2 = 0.5,1,2; which displays 
the strange of interaction between matter and dark en- 
ergy. Fig. 2. shows that q has negative values that it 
means the universe is in an acceleration phase and get- 
ting larger. Since there is a positive value for Hub- 
ble parameter, f2 A increases too. By passing time and 
increasing dark energy density, q approaches to zero, 
therefore the universe moves from acceleration phase to 
deceleration phase. However since we are in expansion 
epoch of universe evolution, f2 A is still getting larger 
which causes q grows rapidly (and as we expect larger 
values of interaction coupling constant makes q to grow 
up faster). Reaching energy density to critical energy 
density, the universe is forced to stop expansion and 
change its evolution direction and start to contract in 
turnaround point. In contraction epoch, the universe is 
getting smaller more and more with a large value of q 
at the beginning of contraction, as it has been shown in 
Fig. 3. The Hubble parameter has negative values and 
A decreases therefore dark energy density is diluted, q 
again takes negative values, however the universe con- 
tinue to contract and decrease more and more, until 
it moves to low energy regime and dark energy dilute 
enough in which matter and radiation become domi- 
nant. 

Up to now, from above discussion in high energy 
regime, we found out that NADE can satisfy a phan- 
tom fluid for cyclic scenario. However this is not at all 
our universe, there is also regime related to the low en- 
ergy, where p <C p c . At this regime we have the usual 
Friedmann equation for universe, namely 



n, 



Pi 

3H 2 ' 



(20) 



where i refer to matter and dark energy. One can 
rearrange the standard Friedmann equation as 1 = 
Q m + 17a, from this relation we find out the maximum 
value of fti is one. Now, NADE is taken as dark energy 
component with interaction with matter, so from con- 
servation equation related to the dark energy and with 
the help of above Friedmann equation, one obtain lu\ 
as 



cj a = -1 



3a n 



-36 2 (l + r), 



(21) 



maximum value of dark energy density parameter, 
namely fl\ is one and also when the scale factor is large 
enough to make the second term on the right hand of 
above relation small in which uja stands in phantom 
range. The sort of dark energy as a component of uni- 
verse fluid has been determined, so we want to be aware 
about the evolution of universe in this stage. 
The deceleration parameter is expressed as 



1 



|((1 +uj m )D, m + (1 + cj a )^a); 



(22) 



so in the earlier time when Q\ converge to zero (and 
Q m converge to one) the deceleration parameter is as 
q w 1 for oj m = -| , and q w \ for oj m = that indicates 
a deceleration universe. Whereas in the late time when 
fl\ approaches to one , the deceleration parameter has 
negative sign that introduces an acceleration universe. 



3 Validity of second law of thermodynamics 

So far, we investigated whether NADE can produce a 
phantom dark energy for cyclic universe. Now, it can be 
valuable to study validity of generalized second law of 
thermodynamic in the model and consider which con- 
dition can satisfy generalized second law. Since loop 
quantum gravity correction enters in Friedmann equa- 
tion which plays an important role in the evolution of 
universe, especially in high energy regime, it is not suit- 
able one straightforward extends the standard treat- 
ment o f thermodynam ic. For this reason we use the re- 



sult of (ILi et al. 112003 ) which thermodynamic has been 



3H 2 



P 



PA- 



investiga ted in loop quan tum cosmology there. The au- 
thors of (|Li et al. 1120091 ) have performed their work by 
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taking a dynamical apparent horizon as universe hori- 
zon. It has been demonstrated that one can have 



dE = TdS + WdV. 



(23) 



However note that, here E = p e ffV = p [1 — j-J V 
which p e ff is effective energy density and W = \{p e ff — 
P eff ) (in which P eff =p(l-^j-^is effective 
pressure). The rela tion can be very useful in this step 
of the work. Like (|Li et al. 1 120091 ) . we take apparent 
horizon 

1 



Ra = 



H' 



From considering the validity of generalized second law, 
we should sum variation of entropy for both of horizon 
and fluid inside the horizon. The horizon entropy 



S h = 2irR A R A 



(24) 



For the relation (|2"5|) . the entropy variation of fluid in- 
side the horizon can be acquired as 



TrSj = Vp 



'eff 



(Peff + Peff) 



V 



(25) 



Note that, there is an equilibrium between the temper- 
ature of horizon and matter inside the horizon. Us- 
ing modified Friedmann equation and conservation re- 
lations and also with the help of the definition of hori- 
zon and energy density parameters, the total entropy 
variation is resulted as 



Ra < 



3(1 + cj a )(1-2«) ^ 

ill 



Ti 



However, if a > ^, the total entropy variation is always 
positive. 

In low energy density, tends to f^A where < Qa < 
1, also a approaches to zero, so one arrives at 



3/3 4 

= 3?T/3 ( R A + — - — 



(28) 



where j3 = ((1 + u m )Q rn + (1 + wa)^a) ■ In radiation 
and matter dominant, where Q\ converges to zero, (3 is 
positive and we should have 



Ti > 



3/3 



Ra 



to generalized second law be valid. In radiation and 
matter dominant uj m is equal to ^ and respectively. 
This expresses that the right term is positive. Passing 
time, the dark energy eventually dominates and fl\ in- 
creases to one, so (3 < 0. The generalized second law 
leads us to the below condition for temperature 



T, < 



4 -3/3 
Ra 



and since /3 is negative we can have a positive value for 
temperature. So the generalized second law of thermo- 
dynamics can be valid, and it can be a confirmation for 
the model. 



3tt((i + u> m )n m + (i + cj A )nty (i - ^ 
+ 1- ((i + cj m )n m + (i + cj A )^ A ) 



X 1 



2p 



(26) 



Now, we turn our attention to the high energy regime 
where Q m can be ignored against to the larg value of 
Q\ . Define a = -2- in which < a < 1, then the total 
entropy variation is reorganized as 



3tt(1 +w a )(1 - 2a)n\ 



x \ Ra + y i {1 + lua){1 ~ 2a)n ^ ~ ^}' (27) 

where (1 + uia) < 0. If a < ^, for making valid the 
second law, there should be 



4 Conclusion 

During this work, in brief, it was shown that both of 
ADE and NADE, without interaction with matter in 
both high and low energy regime, can not satisfy phan- 
tom dark energy which it is a necessary condition for 
cyclic universe scenario. So, as it was the main goal of 
the work, an interaction between the fluid components 
of the universe was supposed. At first step, we turned 
our attention to the high energy regime, and we realized 
that , in ADE case for providing a phantom fluid there 
must be a large value for coupling constant b which it is 
not compatible with the results of previous works, but 
NADE can produce a phantom fluid for cyclic universe 
while dark energy density parameter, namely ilJv, be- 
haves as where £ < 2. It was explained that, in 
expansion phase, equation of state parameter of dark 
energy increases by decreasing r = However in 
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contraction phase, the quantity r increases, so ui\ de- 
creases. In contraction phase with attention to this fact 
that dark energy stands in phantom range, smaller val- 
ues of wa express that dark energy density decreases 
faster and let universe to enter matter dominant phase. 
Moreover larger value of interaction coupling constant 
produce smaller valus for ui^. 

Considering differential equation of fi A exhibited that 
after passing the universe from turnaround point (in 
contracting phase), fi A decreases and after a while the 
universe is putted in low energy regime, then matter 
and radiation are be allowed to be dominate. By com- 
puting deceleration parameter, it was found out that 
there is negative values for q which shows that the uni- 
verse is in an accelerated expansion phase. Increasing 
dark energy density enters q in positive area and it is 
made larger. Then, the universe moves from an ac- 
celerated expansion phase to a deceleration phase and 
it is slowing down until it stops expansion and starts 
contraction in turnaround point. So it is realized that 
the model is able to create a deceleration phase and 
start contraction after an accelerated expansion phase. 
In next step, we investigated the low energy regime of 
universe where it can be possible to ignore loop quan- 
tum correction in Friedmann relation and go back to 
usual form. Investigation shew us that, in low energy 
regime, lu\ stands in phantom range. Also, considering 
deceleration parameter displayed that universe expan- 
sion decelerates in both of matter and radiation domi- 
nant, where J1a is small enough, and accelerates by in- 
creasing J7a toward one. In the last part of the work, we 
studied the validity of generalized second law of ther- 
modynamic in both case of high and low energy. In 
high energy it was determined that when a = y- > | 
the total entropy variation is always positive. At low 
energy, we shew that to make valid the law, a bound 
for temperature is imposed. 
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